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ABSTRACT 

The power spectrum of the two-degree Field Galaxy Redshift Survey (2dFGRS) 
sample is estimated with the discrete wavelet transform (DWT) method. The DWT 
power spectra within 0.04 < k < 2.3hMpc~ 1 are measured for three volume-limited 
samples defined in connective absolute magnitude bins —19 ~ —18, —20 ~ —19 and 
—21 ~ —20. We show that the DWT power spectrum can effectively distinguish ACDM 
models of as = 0.84 and as = 0.74. We adopt maximum likelihood method to perform 
three-parameter fitting with bias parameter 6, pairwise velocity dispersion a pv and 
redshift distortion parameter (3 = Q m 6 /b to the measured DWT power spectrum. 
Fitting results denotes that in a as — 0.84 universe the best fitted Q m given by the 
three samples are consistent in the range 0.28 ~ 0.36, and the best fitted a pu are 
398^27, 475^29 an d 550 ± 20km/s for the three samples, respectively. However in the 
model of as = 0.74, our three samples give very different values of VL m . We repeat the 
fitting by using empirical formula of redshift distortion. The result of the model of 
low as is still poor, especially, one of the best value a pv is as large as 10 3 km/s. The 
power spectrum of 2dFGRS seems in disfavor of models with low amplitude of density 
fluctuations. 
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1 INTRODUCTION 

The present clumpy structures indicated by galaxies on large 
scales are evolved from the very small density fluctuations 
in the early era of the universe. The amplitude of the fluctu- 
ation is fundamental to understand the structure formation. 
A remarkable success of modern cosmology is that the am- 
plitude of mass fluctuations detected by the anisotropy of 
cosmic microwave background radiation is in excellent agree- 
ment with the analysis of the galaxy clustering at low red- 
shifts. Recently released WMAP third year data (WMAP3) 
refines most results of cosmological parameters given by 
the WMAP 1st year data. However, the fluctuation ampli- 
tude smoothed in a spherical top hat windo w of radius of 
8h _1 Mpc is found as small as <rg = 0.74^06 iSpergel et alJ 



2006), which is significantly lower than erg = 0.84 ± 0.04 
of the WMAP 1st year data. The new result of erg is 
a challenge to the cosmological parameter determinations 
from samples of galaxies and galaxy clusters, most of which 
yield as — 0.9 — 1 if the matter content of the universe 
fi m < 0.3 (e.g. iReiprich fc Bohringerll2002t IrToekstra et alJ 



2001 iRefreeier et al.| 



ifregicr et al.1 120021: IVan Waerbeke et al.l 
alJl2003l:lBaficall fc Boddl2003HSeliak et alJ 



Bacon et 



2002 



2005 



Viel fc Haeh nelt 2006). 



The problem motivates us to revisit the constrains on 
as from the power spectrum of the sample of the two- 
degree Field Galaxy Redshift Survey (2dFGRS). The fi- 
nal released spectroscopic catalog of the 2dFGRS contains 
221414 galaxies with good redshift quality Q > 3 and cov- 
ers approximately 1800 square degrees of the sky. It is a 
good sample for studying the fluctuations of cosmic mass 
field on large scales in linear regime as well as on scales 
in nonlinear range. Moreover, the 2dFGRS team has made 
detailed analysis of the Fourier power spectrum, the two 
point correlatio n functions and relevant cosmological pa- 
rameter fitting iPercival et alJ l200ll: iNorberg et alJ feoOll: 



Peacock et alJ 120011: iNorberg et alJ l200lt lHawkins et all 
20031: iPercival et alJ2004:ICole et alJ2005lh Thev found that 
Q m ~ 0-3 or less JPeacock et all200lHCole et al.ll2005lh and 

the best fitting value of as is ~ 0.95 if one takes Q m ~ 0.3 
JPercival et alJl2004) . which is substantially different from 
WMAP 3. 
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In the linear regime ag only controls the overall am- 
plitude of the power spectrum, and it is degenerated with 
the linear bias parameter b. Power spectrum on small scales 
is more effective to constrain as than on large scales since 
the nonlinearity of the power spectrum is directly reflected 
by the value of erg. Most measurements of the 2dFGRS 
power spectrum are on scales of k < 1/iMpc -1 while cos- 
mological parame ter estimation is performed on scale of 
k <~ 0.2/iMpc" 1 JPeacock et alJl20fm IPercival et alll20mt 



iTeemark etail 120021 : ICole et alJ 120051) . With the estima- 
tor based on the dis c rete w avelet transformation (DWT) 
flYang et al.ll200llalbL 120021) . we are to analyze the power 
spectrum of the 2dFGRS sample on scales down to k ~ 
2/iMpc -1 . In this scale range, the redshift distortion of DWT 
diagonal mode power spectrum can be easily approximated 
and the al iasing effect is exa ctly eliminated by the DWT 
algorithm <Fang fe Fendl200fJl . 

The paper is organized as follows. In Section 2, the 
DWT power spectrum estimator is briefly introduced. Sec- 
tion 3 describes the construction of samples. The Section 
4 demonstrates robustness and accuracy tests on the DWT 
power spectrum estimator. Section 5 lists our fitting results. 
Conclusions are stated in Section 6. 



2 DWT POWER SPECTRUM ESTIMATOR 

The method of measuring galaxy power spectrum with the 
multi-resolution analysis of discrete wavelet transformation 
has been developed in the last decade (e.g. U ^andofcRmal 
199.4 Il99d: iFang fc Fenl 1200(1 lYang et alJ l2001bl. 12002 
Zhan fc FanJ2003T ) . A brief summary of the method is given 
here, more details are written in the Appendix. 



Fang (2006, hereafter paper I) that the power of the fluctu- 
ations on the modes with the scale index j can be estimated 
by 
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The term I? is the mean power of j modes measured from 
the observed realization n 9 (x), and Nj is the power on j 
modes due to the Poisson noise. For a volume-limited survey, 
the mean galaxy density n g is independent of the redshift. 
The Poisson noise power is thus simply l/n 9 . Pj is usually 
referred to as the DWT power spectrum. 

The DWT power spectrum Pj is related to Fourier spec- 
trum P{n\, 712, 713) by 
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|V(ni/2 J1 )V(n2/2 J2 )V(n 3 /2 J3 )|' ! P(n 1 ,n2,n3) 



where ij)(n) is the Fourier transform of the basic wavelet 
tp(x). Since |-i/j(n)| 2 is a high pass filter in the wavenum- 
ber space, Pj is banded Fourier power spectrum. If the cos- 
mic density field is isotropic, the Fourier power spectrum 
P(ni, ri2,nz) depends only on n — \Jn\ +n| + n%. Eq.(6) 
are exact for homogeneously random fields, either Gaussian 
or non-Gaussian. 



2.1 DWT Power Spectrum 

The observed galaxy number density distribution is 



,(x) = ^2 WmS D (x - X„ 



(1) 



where A^ 3 is the total number of galaxies, x = (xi, £2, £3) 
is a 3-dimensional position vector, x m is the position of the 
m th galaxy, w m is its weight, and 8 D is the 3-D Dirac 5 
function. For an observed sample, n 9 (x) can be regarded 
as a realization of a Poisson point process with intensity 
n(x) = n 9 (x)[l + <5(x)], where n 9 (x) is selection function, 
and <5(x) is the density contrast. 

In terms of the DWT decomposition, the galaxy field is 
described equivalently by variables defined as 



(2) 



= / <5 g (x)^j,i(x)dx , 



where <5 9 (x) = [n 9 (x)/n 9 (x)] — 1, and ?/jj,i(x) is the basis of 
the DWT decomposition, where index j = (ji, 32, js) stands 
for the scale, and 1 = {h,h,h) for position (see Appendix 
A). Since the orthogonal-normal bases i/'j i i(x) are complete, 
all second order statistical behavior of the field can be de- 
scribed by The goal of power spectrum measure- 
ment is to estimate the power spectrum of the density fluctu- 
ations <5(x) = [n(x)/n(x)] — 1 from the observed realization 
8 g (x) = [n 9 (x)/n 9 (x)] — 1. It has been shown in Fang & 



2.2 DWT algorithm of redshift-distortion 

The DWT power spectrum depends on the scale and shape 
of the DWT mode ^^(x), it is sensitive to distortion of the 
shape of the field. It is necessary to establish the mapping 
from the redshift space to the real space. The mapping is 
attributed to bulk velocity and pairwise peculiar velocity. In 
the linear treatment of bulk velocity, the redshift-distorted 
DWT power spe ctrum, is relat ed to the real space power 
spectrum Pj by jYang et al.ll2002h 

P j s = & 2 (l+/3S i ) 2 Sf v P J , (7) 

in which f3 = f2^ 6 /6, b is the linear bias parameter. Sj of 
Eq.(7) is the linear redshift distortion factor. For a cubic 
box of Li = L2 = Lz = L, 



2J1+J2+J3 



E 



ni + ni + ni 



(8) 



|V(ni/2 J1 )^(n 2 /2 J2 )^(n 3 /2 M )| 2 



For diagonal modes ji = J2 = J3 = j, Sj,jj = |. The factor 
Sf v in Eq.(7) is the pairwise velocity dispersion factor. In 
the plane-parallel approximation, if the direction j'3 is chosen 
to be the line of sight, we have Sf v = [s?"] 2 , with 

L f |v3(na/2-)j 2 exp[-4 i (^) 2 ]- 0) 



2J3 
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Obviously, the factor (l+/3Sj) 2 corresponds to the linear 
redshift distortion, and Sf v is the nonlinear redshift distor- 
tion caused by the pairwise velocity dispersion. Although 
Eq.(7) is given by the linear approximation of bulk velocity, 
N-body simulation points out that the mapping of Eq.(7) 
works well till scale k ~ 2/iMpc -1 , also because Pj is weakly 
affected by non-linear clumps of the density field. In general 
for non-volume limited samples, selection function shall be 
taken into account to model redshi ft distortion which brings 
in high order correction to Eq.(7) iYane et alJl20o3) . 



3 SAMPLE CONSTRUCTION 



Samples used in our anal ysis are constru c ted b asically in 
the same way as that in iPan fc SzapudU (120051) . We cre- 
ate volume limited samples from th e 2dFGRS spectro - 
scopic catalog of the final data release llColless et al] l2003). 
which c ontains 221414 gala xies with good redshift quality 
Q > 3 iColless et alJl200ll) . We exclude the ancillary ran- 
dom fields, leaving the two major contiguous trunks: one 
near the South Galactic Pole (SGP) covering approximately 
-37°.5 < S < -22°.5, 21 h 40 m < a < 3 h 40 m and the other 
around the North Galactic Pole (NGP) defined roughly by 
-7°.5 < S < 2°5, 9 h 50 m < a < 14 h 50 m . In order to get 
maximum number of galaxies while keep an uniform sam- 
pling rate to guarantee the fairness of our statistics, we try 
different values of completeness / (/ is defined as the ratio 
of the number of galaxies with redshift to the total number 
of galaxies contained in the parent catalog): fields with com- 
pleteness less than the chosen value is cut off, and fields with 
higher completeness is diluted to match the sampling rate. 
We find / = 0.738 is the optimal value. The final parent sam- 
ples are thus restricted by completeness / > 0.738, and ap- 
parent magnitudes limits in photometric bj band with bright 
cut of mi, j — 15 and faint cut of median value of ~ 19.3 wit h 
some small variation specified by masks llColless et al .120031) . 

Volume limited sub-samples are built from the parent 
sample by selecting galaxies in specified absolute magnitude 
ranges. A bsolute magnitudes ar e calculated with k + e cor- 
rection in lNorberg et all i2002f) . Our analysis focus on the 
three sub-samples defined in absolute magnitude Mt r bins of 

-19 ~ -18, -20 ' 19 and -21 r 20. Basic parameters 

of the three volume limited samples of 2dFGRS are summa- 
rized in Table in which lists range of redshifts z m i n -z max , 
range of comoving distances dmin- d max , numbers of SGP 
and NGP galaxies, and the mean densities n. Comoving dis- 
tances are calculated from redshifts z in the ACDM universe 
with J1a = 0.7 and Q m = 0.3. 



4 NUMERICAL TESTS OF THE DWT POWER 
SPECTRUM ESTIMATOR 

In this section, we will test the DWT power spectrum esti- 
mator with Poisson samples and samples from N-body sim- 
ulation. Nine realizations of Poisson samples are produced 
in box with side L = 239.5 h _1 Mpc with 256 3 points. All 
the measured DWT power spectra of these samples are plot- 
ted in the top panel of Figure Q To test the stability of the 
DWT estimator we calculate the diagonal DWT power spec- 
trum, Pj,j,j for each realization, and then computing their 



or 




J 

Figure 1. Top: Diagonal DWT power spectra of nine realizations 
of Poisson samples. Each sample is generated in a box with size 
of 239.5 h — 1 Mpc and 256 3 particles. The relation between j and 
physical scale is I = 239.5/2 J h — 1 Mpc. Bottom: Ratios of diag- 
onal DWT power spectra of nine realizations of Poisson samples 
to their mean spectrum. 



mean Pj,j,j. Pj.j.j / Pj,j,j are shown in the bottom panel 
of Figure We can see that at j < 2(i.e. scales larger 
than 119.75/i _1 Mpc), there are as large as 50% variances 
in the diagonal DWT power spectrum. Thus we will not 
use data points with j < 2. At small scales, or large j, the 
DWT estimator gives reliable results. This is because the 
aliasing effect is effe ctively suppressed in the DWT analysis 
jFang fe Fendl2000l) . 

Then, in order to test the geometric effect of samples on 
estimation of the power spectrum, we cut one of the Poisson 
sample into three sheet-like sub-samples as 60.0 x 239.5 x 
239.5, 20.0x239.5x239.5, and 20.0x60.0x239.5 (h _1 Mpc) 3 . 
Furthermore, a fourth sub-sample is constructed from the 
20.0 x 239.5 x 239.5 (h^Mpc) 3 sub-sample by cutting off 
three parallel cylinders with radius of 5.0, 10.0 and 20.0 
/i _1 Mpc respectively. In Figure |5J we plot the ratios of the 
DWT power spectra of each sub-sample to their parent ran- 
dom sample P? j j/ Pj,j,j- For j = 2 and 3 (or on the scale 
of 239.5/2 h _1 Mpc and 239. 5/2 2 h _1 Mpc), the scatters in 
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Table 1. Volume limited sub-samples of 2dFGRS. 
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Figure 2. The ratios of DWT power spectra of cutted Poisson 
sub-samples to the parent Poisson sample. The parent Poisson 
sample is generated in a box with size of 239.5 h _1 Mpc, with 
256 3 points. The sub-samples are defined in box of 60.0 X 239.5 X 
239.5/i~ 1 Mpc (dash-dot-dot line); 20.0 X 239.5 X 239.5ft.- 1 Mpc 
(dash line); 20.0 X 60.0 X 239.5h~ 1 Mpc (dash-dot line); and 20.0 X 
60.0 X 239.5/i — x Mpc with three cylinders chopped (dot line). 
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Figure 3. DWT power spectrum of Virgo simulation sample. 
The dash line on the right is measured with the box size of 479.0 
h~ 1 Mpc. Shifting the right dash line by one unit along j-axis, we 
get the dash line on the left. The dot line is the power spectrum 
measured in box of size 239.5 /i _1 Mpc. The solid line is th e non- 
linear power spectrum from the formula of ISmith et alJ 120 
with C m = 0.3, n A = 0.7 and erg = 0.9. 



those spectra can be larger than 50%, and the ratios are 
randomly distributed with variance of the order of unity. 
Actually, such significant scatters result from small number 
of modes on j < 3. For j > 3, differences between those 
spectra are negligible, the DWT power spectrum estimator 
is well independent of sample geometry on small scales. 

To further test the reliability of the DWT power spec- 
trum estimator, we measured the DWT power spectrum of 
the Virgo simulation which is of ACDM cosmology with 256 3 
particles in box of size 239.5/i _1 Mpc. Figure compares the 
theoretical DWT power spectrum with that estimated from 
the simulation. The theoretical DWT power spectrum is cal- 
culated from Eq.(6) with n onlinear power spect rum from the 
accurate fitting formula of lSmith et alJ i2003l) . Clearly, the 
theoretical spectrum and measurements are in good agree- 
ment on scales less than 239. 5/2 4 ~ 15/i~ 1 Mpc. The test 
shows that the DWT power spectrum estimator can per- 
fectly recover the original power spectrum on small scales. 

The final test is to measure the DWT spectrum of the 
Virgo sample in a cubic window of side 479.0/i _1 Mpc which 
is twice of the simulation box size. Theoretically the power 
at scale j in the box of side 479.0/i~ 1 Mpc corresponds to 
that at scale j — 1 in the box of side 239.5/i _1 Mpc. It is 
clearly seen in Figure |3] that the spectrum measured in the 
479.0/i _1 Mpc box at j exactly equals to that measured in the 
239.5/i _1 Mpc at j - 1. The DWT estimator is independent 



of the size of the window box. One can choose the size of box 
freely and put the sample in the box wherever one like. Note 
that the difference of spectra at j — 2, 3 are greater than 
other data points. Therefore, in the following analysis of the 
2dFGRS catalogs, the first two data points in our spectra 
are discarded. 



5 DWT POWER SPECTRUM OF 2DFGRS 
SAMPLES 

5.1 The diagonal DWT power spectrum 

To achieve the largest possible volume, the NGP and the 
SGP regions are measured together in a window box of size 
1403.0ft _1 Mpc. The size of the window is sufficient to cover 
all three volume limited sub-samples. The filling factor of 
each sub-sample is 0.07%, 0.24%, and 0.68% respectively. 

To estimate error bars of the DWT power spectra, we 
created mock volume limited samples from the 22 mock 
galaxy catalogs which are extracted from the Hubble vol- 
ume s imulation 1 . Details of mock catalogs are in lCole et alJ 
( 1998). The set of the mock samples used here is the Lamb- 
daCDM04. The 22 mock samples are filtered with the same 

1 http:/ /star-www. dur.ac.uk/~cole/mocks/hubble. html 
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Figure 4. DWT power spectrum of 2dFGRS volume-limited 
samples and theoretical spectra. The size of the cubic box is 
1403.0 /i — 1 Mpc. The relation between j and physical scale is 
I = 1403.0/2^ h- 1 Mpc. The power at large scales are suppressed 
by random motion of galaxies. Due to cosmic variance, the error 
bars at j = 3 is large. At j = 9, Poisson noise leads to large error 
bars. 

selection criteria and masks as the real galaxy volume lim- 
ited sub-samples. Error bars of the DWT power spectra of 
the three 2dFGRS sub-samples are approximated by the 1-a 
scattering of their mock samples. 

Short noise is not directly calculated with Eq.(5). In- 
stead, we produce a number of Poisson samples with the 
same geometry as our samples. Numbers of points in these 
Poisson samples are the same as the number of galaxies in 
the sample under analysis. 22 Poisson samples for each sub- 
sample allow us to estimate error bars for short noise sub- 
traction. These error bars are incorporated into the final 
results of galaxy power spectra. 

In Figure 2] we present the measured diagonal DWT 
power spectra of the three volume limited samples, together 
with two theoretical nonlinear spectra of fiat ACDM model 
with parameters A.) Q m — 0.3, = 0.7 and erg = 0.84 
(Model A), and B.) n m = 0.24, Q A = 0.76 and erg = 0.74 
(Model B). The scale range is 0.035 < k < 2.2/iMpc" 1 . 

5.2 The fitting of redshift distorted power 
spectrum 

The differences between the power spectra of model predic- 
tions and the real data shown in Figure 2] are mainly due 
to redshift distortion and bias. We adopt Eq. (7) to fit the 
power spectrum of 2dFGRS samples with the nonlinear real 
sp ace DWT power sp ectrum from Eq. (6) and the formula 
of ISmith et all J2003D . We firstly use the 22 mock samples 
to estimate the correlation between the powers of different 
DWT modes. We construct a covariance matrix 

d n' = J2 Ad ' A 4 > ( 10 ) 

where N Bim — 22 and Adj = di—{dj). The vector d 1 consists 
of elements = Pjjj with j — 3, 9, Pjjj is the power 



spectrum from the 7-th simulation, and (dj) is the mean. We 
found that off-diagonal elements of the covariance matrices 
are always one order of magnitude smaller than diagonal 
elements. Actually it is a typical feature in the DWT anal- 
ysis. The correlations between different modes are highly 
suppressed regardless the field is Gaussian or non-Gaussian 
(Feng & Fang 2005). The quasi-diagonalizing of the correla- 
tion matrix in the DWT decomposition has been extensively 
used for data compression. By virtue of this property, we can 
compute x 2 with diagonal elements only. Namely, we use the 
Chi-Square as our maximum likelihood estimator 

2 ^ [P ] ~pS{a 1 ,a 2 ,az...)] 2 

x =2^ ~2 ■ ( U ) 

in which Pj is observation data and Pf(ai, ai, 03...) is 
the redshift-distorted power of model with parameters 
(ai, <J2, H3-..). We take the Reduced- Chi- Square, which is de- 
fined as Xd.o.f = X 2 K N - M ) where (N - M) is the degree 
of freedom, as our final results shown in tables. 

We aim at detecting the influence of as on the power 
spectrum. Two fiducial ACDM models are considered here: 
A. erg = 0.84, and B. <rg = 0.74. We select the linear bias 
parameters b, redshift distortion parameter (3 (or fi m ) and 
pairwise velocity variance a pv as fitting parameters. Other 
parameters of model A and B are the same. The parame- 
ter space (b, 13, a pv ) is divided into a 20 x 20 x 20 grid. The 
first run of fitting is performed on very crude grids in broad 
parameter space to locate the region of best \ 2 ■ Then we de- 
crease the volume of parameter space centered in this region 
with finer grid to obtain the three dimensional probability 
distribution functions (PDF) of (b,/3,a pv ). After integrating 
over two of the three axes in the parameter space, we have 
the marginalized PDF for each parameter. 



5.3 Fitting results 

The best fitted power spectra of both models are very close 
to each other, so only the fitted power spectra of model A 
are demonstrated in Figure The estimated parameters of 
the model A and B are tabulated in Tables [5] and [3] respec- 
tively. As an example, marginalized PDFs of the parameters 
b, f3, <j pv for the sub-sample —21 ~ —20 are in Figure|S| Fig- 
urc[7|gives the PDFs of Q. m for all the three sub-samples. 

In Table [5] we can see that three sub-samples offer 
about the same estimation of Q m — 0.3, and the b increases 
from 0.75 to ~ 1 with the luminosity of gala xies, which 
is consistent with an alysis of others iNorberg et al.l l200ll : 
iPan fc Szapudll2005l) . The pairwise velocities obtained in 
the three samples are in the reasonable range, and also in- 
crease with galaxy luminosity. It appears that our analysis 
with model A (erg = 0.84) is basically in good agreement 
with previous works, and more importantly, that the DWT 
proves itself an effective tool for parameter estimation from 
galaxy samples. 

Fitting to model B provides very different estimation of 
parameters. As seen in Table 01 and Figure [7] values of f2 m 
given by the three sub-samples are quite different from each 
other. The best values of O m for the sub-samples of — 19 ~ 
— 18 and —20 ~ — 19 are significantly larger than 0.3, which 
is in disagreement with most current measurements at least 
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Table 2. Parameters estimated by fitting 2dFGRS DWT power spectra with Eq. (7) to ACDM model of cr 8 = 0.84. 

M b -51og 10 h & /Cl m b o pv (kms^ 1 ) x' 2 



-J.!> -in 0.66±g"?g /0.31±g"?| 0.75+. , ,! 39S±1? 8.87 

,+0.07 /n o fi +0:i2 noR+OiOS ^r+37 



-20 19 0.621^ /0.36I(, 14 O.86T00? 475+29 16 - 19 

-21 --20 0A3±™t /0.28tg-;g l.OTtgft 550+J° 4.06 



Table 3. Parameters estimated by fitting 2dFGRS DWT power spectra with Eq. (7) to ACDM model of a s = 0.74. 

M b j - 5 log 10 h P /n m b ct p „ (kins" 1 ) \ 2 

-19 --18 0.76l°;<£ /0.4ll°;i° 0.77±g;g| 415+*° 12.24 

-20 — -19 0-73l ( o;lo /0.46+.°; 2 ° - 86 -o:o7 492+.^ 17.99 

-21 — -20 0.35l°;° 2 /0.28l°;o2 i-^toio! 600^^ 4.83 



Table 4. Parameters estimated by fitting 2dFGRS DWT power spectra with Eq. (7) to the model of f2 m = 0.24 and Ha = 0.76. 
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Figure 6. Marginalized distribution of parameters b, (3, and cr pv in fitting the DWT power spectrum of the 2dFGRS sample of M bj £ 
(-20 21) with Eq. (7). 
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Figure 7. Marginalized distributions of Q m in fitting the DWT power spectra with Eq. (7) for 2dFGRS samples with M 6/ in — 18 ~ 19, 
-19 ~ -20 and -20 21. 
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Mb A 



Kh" 1 Mpc) 
43.8 21.9 




a pv : the PDF of model B is highly skewed and very broad 
while the PDF of model A is close to Gaussian. 

The result of model A shown in Table 03 is roughly the 

same as that in Tabled The sub-sample of -19 ' 20 which 

gives a large Q m but still agree with others within error bars. 
The results of model B shown in Table |5] are similar to those 
in Table |3] Q m from all sub-samples are larger than 0.3, es- 
pecially the one from the sub-sample of —19 ~ —20 which is 
unusually as large as 0.731q'27- Meanwhile, the pairwise ve- 
locity variance estimated from the sub-sample of - 



has an extraordinary value of 1014.28l^3 4 5 7 kms 



21 ~ -20 
_1 . These 

impose questions on the prior of as = 0.74. 

Finally, we take ag , b and a pv as fitting parameters, and 
fix Q m = 0.24. The best fitting parameters are in Table |7| 
Again we have the similar results as those listed in Table [I] 
All the values of as are always > 0.8, which suggests that 
the low as is unlikely to match with the 2dFGRS. 



Figure 5. The DWT power spectra in redshift space of the 
ACDM model with trg = 0.84. Parameters b, ft, and fjpy arc 
given by our best fitting to the data points shown in the Fig- 
ure. The three panels are for sample within absolute magnitudes 
-20 ~ -21 (top); -19 ~ -20 (middle) and -18 ~ -19 (bottom). 



at level (e.g. iPeacock et aljEoOll iTeemark et al]l2004h 
Only the sub-sample of —21 ~ —20 yields the Q m « 0.3. 

In order to place constrains from the 2dFGRS on as, 
we repeat the fitting procedure with as, b and a pv as fitting 
parameters, and fi m = 0.24 as prior. Results are shown in 
Table 0] It is very clear that, except the —19 ~ —18 sub- 
sample, the other two sub-samples both give as > 0.9. An 
universe of as = 0.74 is not preferred by the DWT power 
spectrum of the 2dFGRS. 



5.4 Fitting with an alternative formula of redshift 
distortion 

Empirically, the Fourier power spectrum P s (k) in redshift 
space is related to that in real space by dPeacock fc Doddsl 
1994) 



P S (k) 



b 2 P(k)G(y,f3) 



in which y 2 — k 2 a pv , and the function G is 

G (y,P )= 4 e jM [ 3f 3 2 + ^y 2 + ^ ] 



cxp(-y 2 
% 4 



■[/3 2 (3 + 2y 2 )+W 



Substituting Eq.(6) into Eq.(12) yields 



j(b,/3,a pv ) 



(12) 



(13) 



(14) 



where y(j) = k 2 (j)a pv and k(j) is given by Eq. (B12). 

Following the same procedure in Section 5.2, we perform 
a three parameter (b,/3,a pv ) fitting with Eq.(14). Fitting 
results of model A and model B are written in Tables and 
^respectively. PDFs of the three parameters b, f3, a pv for the 
sub-sample of —21 ~ —20 are plotted in Figure 00 and the 
PDFs of fi m for all the three sub-samples are in El Figure Q3 
shows that the models A and B have very different PDF of 



6 CONCLUSION AND DISCUSSION 

DWT power spectra of 2dFGRS samples are measured on 
scales equivalent to 0.04 < k < 2.3/tMpc -1 . We show that 
these power spectra are efficient to test ACDM models with 
high and low amplitudes of mass density fluctuations. The 
model with as = 0.84 finds good support from the 2dF- 
GRS sample, all the best fitting parameters, b, (3, a pv are 
consistent with other works on 2dFGRS. Especially, three 
volume- limited samples gives different b and /3, but approxi- 
mately the same values of fi m = (6/3) 1/0 fi = 0.28 - 0.36. On 
the other hand, the model with as = 0.74 cannot give such 
consistent fitting result, the best fitted fl m from the three 
volume- limited samples are significantly different, and in de- 
viation from ~ 0.3 at least at la level. Moreover, the fitting 
results of a pv are generically large, even reach 10 3 km s _1 . 
Our studies suggest that the power spectrum of 2dFGRS 
disfavors models with low amplitude of mass fluctuations, 
as = 0.74, if other cosmological parameters are given by the 
WMAP3. 

It is found that a pv increases with lumi nosity, which is 
basica lly consistent with the observation of Ijing fc Bornerl 
(2004) though our estimations are lower than theirs. 

Another parameter that will affect the shape of the 
DWT power spectrum is the slope of the primordial fluctu- 
ation spectrum n s . We have been using the scale-invariant 
spectrum, or the Zeldovich spectrum with n s — 1. However, 
we notice that WMAP3 gives n s = 0.95±o:oi9, which is 5% 
smaller. In order to check the influence of a lower value of 
n s on our fitting results, we repeat our fitting for the sub- 
sample -20 < 21 with n s = 0.95. We find b = 1.07±° ° {, 

(5 = 0.44j;°;° 2 (On = 0.29lg;o2), and a pv = 550±?,okm a" 1 . 
Clearly, the change of n s from 1 to 0.95 results in little mod- 
ification to b and a pv , only slightly alters (less than 5%) the 
value of P, or fl m . 

In this paper, we used only the diagonal modes in term 
of vector j, Pj = (ej,i£j,i), andj = Even in the second 

order statistics of the DWT, {ej,i?j,i'), we can have the power 
of off-diagonal modes P- } — (ej.iej.i), and j = {31,32,33)- In 
addition, we have the correlation between modes of (j, 1) and 
(j,l') (1 7^ 1')- It nas been shown that different parts of the 
second order statistics of DWT contains d ifferent informa- 
tion of the random field jYang et alJl200lllah . Possible con- 
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Figure 8. Marginalized distribution of parameters /3, 6, <j pv for the —20 ~ —21 sub-sample, with the empirical rcdshift distortion mapping 
formula of Eq. (12). Fitting with erg = 0.74 give larger b and a pv than with erg = 0.84. 
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Figure 9. Marginalized distribution of S7 m in three parameter fitting with Eq. (12). For all sub-samples, fitting with lower ag gives 
larger value of S7 m . 



Table 5. Best parameters from fitting DWT power spectra to ACDM model with ag = 0.84 and Eq. (12). 



M 6; -51og 10 /i j3 /n m b a pv (kms x ) \ 2 



-19 — -18 0.771°,-°! /°- 34 20 °- 69 -oo9 3151°° 10.91 

-20 — -19 0.73+°;^ /0.56+°;^ 0.97^°;°* 578+^°, 8.46 

-21 20 0.4ll°'o3 /0-28toH 1.15l°'o3 54S -26 4 ' 96 



Table 6. Best parameters from fitting DWT power spectra to ACDM model with ag = 0.74 and Eq. (12). 



M b j - 5 log 10 h f} /H m b a pv (kms 1 ) x' 2 

-19 18 0.76±g-57 /0.38l°-2o 0.74toJ? 336±H 10.61 

-20 --19 0.73l°^ /0.73l°i? 1.13± ,: ° 1 ? 682+« 8.26 

-21 --20 0.40±g;gl /0.39±g;^ 1.42±°;g| 1014l^ 3 4.89 



Table 7. Best parameters erg, b, a pv from fitting DWT power spectra to ACDM model with Q m = 0.24, Qjy = 0.76, and Eq. (12). 



M bj - 


5 log 10 h 


C8 


b 


a pv kms 


x 2 


-19 - 
-20 - 
-21 - 


-18 
-19 
-20 


0-80±g;SS 

0.93+°,;J° 


73+ - 21 
91+°- 15 

97+0-05 
U ' M ' -0.04 


325±f° 
470+" 
412+ 50 


17.00 
15.62 
7.16 
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strains on parameters given by various parts of the second 
order DWT statistics deserve further study. 
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APPENDIX A: DWT DECOMPOSITION OF 
RANDOM FIELD 

For the details of the mathem ati cal properties o f the D WT 



latical properties ot the JJ W 1 , 
lease refer to Mallat i I ?)s< )ailb|) : iMeveil dl992lk iDaubechies 



( 1992), and for physical applications, refer to lEme^^Tliew; 
(1998). For this application, the most important properties 
are 1.) orthogonality, 2.) completeness, and 3.) locality in 
both scale (r) and physical position (a;). Wavelets with com- 
pactly supported basis are an excellent means to analyze 
random fields. Among the compactly supported orthogonal 
wavelets, the Daubechies family of wavelets are easy to im- 
plement. 

To simplify the notation, we consider a 1-D field p(x) 
on spatial range L. It is straightforward to generalize to 3-D 
fields. In DWT analysis, the space L is chopped into 2 J seg- 
ments labeled by I = 0, 1, ...2 3 — 1. Each of the segments has 
size L/2 J . The index j is a positive integer which represents 
scale L/2 3 . The index I gives position and corresponds to 
spatial range lL/2 j < x < (I + l)L/2 j . 

DWT analysis uses two functions, the scaling functions 
<t>j,i(x) = (2 3 /L) 1/2 4>(2 3 /L - I), and wavelets ipj,i(x) = 
{2? I V) 1 ^{2? j L — I). The scaling functions and wavelets 
are given, respectively, by a translation and dilation of the 
basic scaling function <f>(rf) and basic wavelet ip(ri) as 



and 



<t>iA x ) 



1/2 



cj>{2 j x/L-l) 



1/2 



ip{2 J x/L-l). 



(Al) 



(A2) 



The scaling functions play the role of window function. They 
are used to calculate the mean field in the segment I. The 
wavelets ipj,i(x) capture the difference between the mean 
fields at space ranges IL/2 3 < x < (I + l/2)L/2 3 and (I + 
l/2)L/2 3 < x < (Z + l)i/2 J . 

The scaling functions and wavelets tpj t i(x) satisfy the 
orthogonal relations 



(j> jt i{x)(j)j > v{x)dx = S u , 



(A3) 
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J 



ip ji i{x)ip J , >l ,{x)dx = 8 j ^8 l>l ,, 



(A4) 



4> j ,i{x)ip j , tV {x)dx = 0, iij'>j. (A5) 



With these properties, a 1-D random field p(x) can be 
decomposed into 

oo 23'' -1 

p(x) = P 3 {x) + J2 tj'.rtj'Ax), (A6) 



where 



23-1 



P 3 ( x ) = Yl e i^oA x )- 



(A7) 



The scaling function coefficient (SFC) ej t i and the wavelet 
function coefficient (WFC), are given by 



/ 



£j,i = / p{x)4>j,i{x)dx, 



and 



(A8) 



(A9) 



respectively. The SFC e^j measure the mean of p(x) in the 
segment I, while the WFC lj t i measures the fluctuations (or 
difference) of field p(x) at I on scale j. 

The first term on the r.h.s. of Eq.(A6), p j (x), is the field 
p(x) smoothed on the scale j, while the second term contains 
all information on scales > j. Because of the orthogonality, 
the decomposition between the scales of < j (first term) and 
> j (second term) in Eq.(A6) is unambiguous. 



power spectrum of perturbations on scale L/n is then de- 
fined as 



P(n) = |e„ 



(B5) 



This is the power spectrum with respect to the Fourier de- 
composition. 

Similarly Parseval's theorem for the DWT is 



fL oo ,23-1 

±r / Hx^dx^jY. 

L J " j=0 1=0 



(B6) 



Thus, the second order statistical behavior of e(x) can be 
described by the \ij,i\ 2 and one can call \ejj\ 2 the DWT 
power spectrum. 

Comparing Eqs.(B4) and (B6), it is clear that 

~t Sf=o 1 I^j,'| 2 i s a measure of the power of the perturba- 
tion on scales from L/2 J to L/2 j+1 . Therefore, the power 
spectrum with respect to the wavelet basis can be defined 



23-1 



(B7) 



Since the DWT bases ipj t i(x) measure the differences be- 
tween the local mean densities at adjoining scales, the mean 
density on length scales larger than the sample size is not 
needed in calculating The spectrum Eq.(B7) will not 
be affected by the infrared (long-wavelength) uncertainty of 
the mean density. 

With Eqs.(B2), (B6), (A2) and (A4) we can find the 
relation between the power spectra of DWT Pj and Fourier 
P(n). It is 



1 00 

P i=V E \i>W2 j )\ 2 P(n), 



(B8) 



APPENDIX B: 1-D DWT POWER SPECTRUM 

The contrast (or perturbation) of the field p(x) is defined by 



5(x) = 



p(x) - p 



(Bl) 



where p is the mean density of the field. The Fourier expan- 
sion of e is 



8{x) = Y e " e 



n— — oc 



with the coefficients given by 



8{x)e- l2m/L dx. 



(B2) 



(B3) 



Parseval's theorem relates the power for a distribution 
to the coefficients of the Fourier expansion. For the contrast 
this yields 



z j o \e(x)\ 2 dx J2 W 



(B4) 



which shows that the perturbations can be decomposed into 
domains, n, by the orthonormal Fourier basis functions. The 



where ip{n) is the Fourier transform of the basic wavelet 
given by 



i>(n) 



[ i>{n)e 
■Jo 



drj. 



(B9) 



Since wavelet is admissible, i.e. J ip(r))dn = 0, we have 
^(0) = 0. ^(n)! 2 is localized in n-space. |V>(n)| 2 has sym- 
metrically distributed peaks with respect to n = 0. The first 
highest peaks are non-zero in two narrow ranges centered at 
n — ±n p with width An p . Besides the first peak, there are 
"side lobes" in |-0(n)| 2 . However, the "side lobes" is small, 
for instance, for the Daubechies 4 wavelet, the area under 
the "side lobes" is not more than 2% of the first peak. There- 
fore, Pj is a good estimation of the band-averaged Fourier 
power spectrum centered at wavenumber 



The band width is 



71 j — f^p 2 



An = 2 J An„. 



In other words, the relation between k and j is 

log k = (log 2) j - log(L/27r) + log n p . 
For the D4 wavelet, logn p = 0.270. 



(BIO) 



(Bll) 



(B12) 
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APPENDIX C: 3-D DWT POWER SPECTRUM 



where the window function Wh is 



For 3-D random field, the DWT decomposition is based 
on the orthogonal and complete set of 3-D wavelet basis 
{^li,i( x )}i which can be constructed by a direct product of 
1-D wavelet basis as 



(CI) 



where ji = 0,1,2.. (i = 1,2,3) and h = 0...2 5 ' 4-1 . Ob- 
viously, the basis Vj,i( x ) is non-zero mainly in a volume 
Li/2 J1 x x Lz/2 i3 , and around the position (xi 



/2 31 ,x 2 = ; 2 L 2 /2 J 



l 3 L 3 /2 J3 ) (Fang fc Thewj 



Similar to Eq.(B7), the power spectrum on scale j = 
(jl, h, 33) is 



231 +32 +33 



2^1-1 2*2-1 233-1 

E E E 

ii=o ; 2 =o i 3 =0 
For 3-D samples, Eq.(B8) is generalized as 

1 



ej,i 



(C2) 



231+32 +33 



E E E 



ni = -oo ti-2 — — 00 713 — — 00 



(C3) 



|V(m/2 ,1 )V(na/2 ,a )V'(n3/2 ,8 )rP(ni > n a ,iis) 



Because the cosmic density field is isotropic, the Fourier 
power spectrum P(ni, 112,713) is dependent only on 



n l + n 2 + n 3- 



(C4) 



Obviously, the DWT power spectrum is invariant with re- 
spect to the cyclic permutation of index as 



Pii, 



(C5) 



Considering Eq.(C3) and Eq.(C4), we can formally de- 
fine a band center wavenumber rij corresponding to the 3-D 
mode j as 



raj = n P V(2 J1 ) 2 + (2^) 2 + (2«)2. 



(C6) 



For an isotropic random field, the Fourier modes with the 
same n are statistically equivalent. However, the DWT 
modes with the same rij [Eq.(C6)] are not statistically equiv- 
alent, because the DWT modes are not rotationally invari- 
ant. A Fourier mode e -i(^/m^i+n 2 x 2 +n 3 x 3 ) cfm be ob _ 

tained by a rotation of mode e -*(?*/m<*i+"v*3+ n **s) ag 
long as n'i + n 2 2 + n 3 2 = n\ + n\ + n\. However, the DWT 
modes don't have the same property. Generally, one cannot 
transform a mode (ji , j 2 , j 3 ) to (j[ , j' 2 , j' 3 ) by a rotation, even 
when rij ~ riji. Because of different configurations between 
them, the condition rijUji generally does not imply 

Pj = Py- (G7) 

This invariance holds only when (31,32,33) is a cyclic per- 
mutation of (3i,32,3z)- 

With this property, one can define two types of the 
DWT power spectra: 1. the diagonal power spectrum given 
by Pj on diagonal modes ji = j 2 = 33 = j, and 2. off- 
diagonal power spectrum given by other modes. 

From Eq.(C3), the diagonal power spectrum Pj = Pjjj 
is related to the Fourier power spectrum by 
00 00 00 
P 3 E E E Wj(ni,n2,n 3 )P(ni,n2,n 3 ), (C8) 



W ] (n 1 ,n 2 ,n 3 ) ¥ -\i,(n 1 /2 j )^(n 2 /2 1 )^(n 3 /2 :i )\ 2 . (C9) 
with the normalization 

Wj(rai,n2,ra3)dnidra2dra3 = 1. (CIO) 



The window function Wj is localized around rai = n 2 = 
n 3 = n p 2 J . Therefore, the diagonal power spectrum Pj is a 
band-average of the isotropic Fourier power spectrum P(n) 
with the central frequency n — y/3n g 2-' . 

There are two types of modes: diagonal mode with 
ji = j 2 = j 3 ; off-diagonal mode for which the three num- 
bers (31,32,33) are not the same. The DWT estimator can 
provide two types of power spectra: 1. diagonal power spec- 
trum given by the powers on diagonal modes, 2. off-diagonal 
power spectrum given by the powers on off-diagonal modes. 
Because the two types of modes have different spatial invari- 
ance, the diagonal and off-diagonal DWT power spectra are 
very flexible to deal with configuration-related problems in 
the power spectrum detection. For off-diagonal modes, one 
can also calculate the linear non-diagonal DWT power spec- 
trum Pj 1 ,j 2 ,33 v i a Eq.(C2). However, in this case, Pj lt j 2 ,j 3 
cannot simply be identified as a band average of the isotropic 
Fourier power spectrum P(n) centered at n = rij. Neverthe- 
less, rij is useful to calibrate the physical scale of a given 
j- 



